Universal Relations for Identical Bosons 
from 3-Body Physics 
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Systems consisting of identical bosons with a large scattering length satisfy universal relations 
determined by 2-body physics that are similar to those for fermions with two spin states. They 
require the momentum distribution to have a large-momentum l/fc* tail and the radio- frequency 
transition rate to have a high-frequency I/lu^^'^ tail, both of which are proportional to the 2-body 
contact. Identical bosons also satisfy additional universal relations that are determined by 3-body 
physics and involve the 3-body contact, which measures the probability of 3 particles being very 
close together. The coefficients of the 3-body contact in the 1/fc^ tail of the momentum distribution 
and in the 1/cj^ tail of the radio- frequency transition rate are log-periodic functions of k and uj that 
depend on the Efimov parameter. 

PACS numbers: 31.15.-p,34.50.-s, 67.85. Lm,03.75.Nt,03.75.Ss 
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Strongly interacting systems present a challenging 
problem in theoretical physics. Some of the simplest 
of such systems consist of particles with short-range in- 
teractions and large scattering lengths. They arise in 
almost all branches of physics, including atomic, con- 
densed matter, high energy, and nuclear physics. Ul- 
tracold trapped atoms with large scattering lengths are 
particularly pristine examples of such systems. In addi- 
tion to the exquisite probes that are available in atomic 
physics, the ability to control the scattering length by 
using Feshbach resonances makes this dimension of the 
system accessible experimentally. 

The simplest many-body systems of particles that 
interact through a large scattering length consist of 
fermions with two spin states. In the past decade, such 
systems have been the subject of intensive investiga- 
tions, both theoretical and experimental, using ultracold 
trapped atoms [l|. A powerful tool for studying these 
systems is universal relations that are determined by 2- 
body physics but hold for any state of the system 0. 
These relations involve the contact, an extensive prop- 
erty of the system that measures the probability for a 
pair of particles in the two spin states to be very close 
together. The first such relations were derived by Shina 
Tan 's'l . Universal relations were subsequently derived for 
radio-frequency spectroscopy photoassociation Q, 

structure factors [8| , and correlation functions related to 
the viscocity An exciting recent development is the 
experimental confirmation of some of the universal rela- 
tions [loj . 

In this letter, we present universal relations for iden- 
tical bosons with large scattering length. Like those for 
fermions with two spin states, these new universal rela- 
tions involve the 2-body contact. They also involve the 
3-body contact, an extensive property of the system that 



measures the probability for triples of identical bosons to 
be very close together. 

For identical bosons in the zero-range limit, there are 
two interaction parameters: the large scattering length a 
and an Efimov parameter that is defined below 11|. 
Observables can depend only log-periodically on with 
discrete scaling factor e'^/*° « 22.7, where sq ~ 1.00624 is 
the solution to a transcendental equation. In the unitary 
limit a — ±oo, there are infinitely many Efimov trimers 
with an accumulation point at the 3-atom threshold. The 
parameter k„ can be defined by the trimer spectrum near 
the threshold in the unitary limit: — (e^^'^/'''')"fi.^K^/m, 
where n is an integer. 

The 2-body contact C2 and the 3-body contact C3 for 
a state with energy E can be defined operationally in 
terms of derivatives of the energy: 
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For a many-body state at non-zero temperature, the 
derivatives should be evaluated at fixed entropy. The 
normalization of C2 in Eq. ((Ta| has been chosen so that 
the tail of the momentum distribution at large wavenum- 
ber k (given below in Eq. ([2])) is C^jk'^- The coefficient 
on the right side of Eq. (ITal) then differs by a factor of 
1/2 from that for fermions with two spin states. The 
normalization of C3 in Eq. (jlbp has been chosen so that 
the 3-body contacts for the Efimov trimers in the uni- 
tary limit are (e~^'^/^'')"K^. Werner and Castin have ex- 
pressed the derivative of E in Eq. ()lbp in terms of the 
3-body Schrodinger wavefunction at small hyperradius 
0. 
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The importance of the contacts C2 and C3 is that there 
are other properties of the system that depend hnearly 
on these quantities with universal coefhcients that are 
determined by few-body physics. One of these universal 
relations gives the tail of the momentum distribution at 
large wavevector k. We normalize the distribution n(k) 
so that / (Pk n{k) / {2tt)^ is the total number of atoms. 
We will show below that the tail of n{k) can be expressed 
as 



n{k) 
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where F{k) is a universal log-periodic function: 
F{k) = A sin[2so ln(fc/K») + 20]. 



(3) 



The numerical constants are A = 89.2626 and (f> — 
—0.669064. The dependence on the state enters only 
through the contacts C2 and C3 . Their coefficients can be 
determined by matching the tail of the momentum distri- 
bution for any convenient state. Werner and Castin have 
calculated C2 and F{k)C^ for an Efimov trimer in the 
unitary limit a — > ±cxd [12]. The 2-body contact for the 
trimer with binding momentum t\ii^ IS C2 = 53.0972 K*. 
By our definition in Eq. (jlb| . the 3-body contact is 
C3 = kJ. We determined A and 4> by matching the pre- 
cise results of Ref. [12] • A connection between the 1/k^ 
tail in Eq. ([2|) and the derivative of the energy in Eq. ()lbp 
was conjectured in Ref. [l^ . 

Another universal relation for identical bosons is an 
Energy Relation that expresses the energy of a system in 
terms of n(fc), C2, and C3. The total energy E is the 
sum of the kinetic energy T, the interaction energy [/, 
and the energy V associated with an external trapping 
potential. The kinetic energy can be expressed as an in- 
tegral over n{k). Because of the k~'^ tail in Eq. ([2]), the 
integral is linearly ultraviolet divergent. If the fc""* tail 
were subtracted, the integral would be ultraviolet finite, 
but it would still depend log-periodically on the ultravio- 
let momentum cutoff because of the k~^ tail. The linear 
and log-periodic dependence on the ultraviolet cutoff are 
both cancelled by U . The Energy Relation, which is de- 
rived below, expresses the sum T + U m. a form that is 
insensitive to the ultraviolet cutoff: 
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where /o = -8.42427. The lower hmit fc > fco in the 
k~^ subtraction term avoids an ambiguity in the value 
of the integral over the infrared region. The dependence 
on the arbitrary wavenumber kg is cancelled by the re- 
maining C3 term, whose coefficient depends explicitly on 



kg. The universal constant /o in Eq. (|4]) was determined 
by matching an expression for the energy of an Efimov 
trimer in the unitary limit derived by Werner and Castin 

Tan's Energy Relation for fermions with two spin 
states is similar to Eq. Q except that there are no terms 
proportional to C3 [3[. Combescot, Alzetto, and Ley- 
ronas proposed that such a relation should also apply to 
identical bosons 13|. Werner and Castin demonstrated 



that such a relation does not hold for an Efimov trimer in 
(2) the unitary limit 12[. Our universal relation in Eq. 



demonstrates that it fails for any state for which the 3- 
body contact C3 is nonzero. 

Another universal relation is the Virial Theorem 
for identical bosons trapped in the harmonic potential 
V{r) — ^muj'^r^, which was derived by Werner [14!]: 



T+U-V=- 



77 C2 — C3. 
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This can be derived by using dimensional analysis, which 
implies that the differential operator 2ujd/dijj — ad/da + 
Kh,i9/9k* is equal to 2 when acting on the total energy 
E = T + U + V . The Feynman-Hellmann theorem and 
the definitions in Eqs. ([T]) imply that the three partial 
derivatives are proportional to V, C2, and C3, respec- 
tively. 

One of the most important probes of ultracold atoms 
is radio-frequency (rf) spectroscopy, in which an rf sig- 
nal is used to transfer atoms to a different spin state. 
In the case of fermions with two spin states, there are 
universal relations that provide sum rules for the rf tran- 
sition rate [4] and control its high-frequency tail d, Q . 
In the case of identical bosons, the universal relations 
for rf spectroscopy also involve the 3-body contact. The 
high-frequency tail of the rf transition rate is 
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where f2 is the Rabi frequency associated with the transi- 
tion. The transition rate is normalized so that it satisfies 
the sum rule Jdw r(w) = ttSI^A^, where N is the number 
of identical bosons. The log-periodic function Grf(w) is 
calculated below: 

G'rt(w) ^ B1+B2 sin[so Inimuj/hnl) + 20i.f], (7) 

where Bi = 9.23, B2 = —13.6, and = 1.33. 

Quantum field theory is a particularly powerful formal- 
ism for deriving universal relations [15l] . The universal 
zero-range limit for identical bosons can be described by 
a quantum field theory with atom field ijj. The Hamil- 
tonian density consists of the kinetic term for ip and the 
interaction term 
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where d = tptl) and t = jpipip are local composite opera- 
tors. We set 7i = 1 for simplicity. To obtain scattering 
length a and Efimov parameter k^, the bare coupling 
constants g2 and 173 must be tuned as functions of the 
ultraviolet momentum cutoff A 16|. If we use a sharp 



cutoff on the momenta of virtual particles, the bare cou- 
pling constants must be chosen to be 
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J 



87r/[l/a-2A/7r], 
~%gl{H + J/aK)/K\ 
h^{C-saS)/{C + s^S), 
[jo+Ji(2^C)+j2(C2- 



(9a) 
(9b) 
(9c) 

S^)]/{C + s^SfiM) 



where C = cos[so ln(A/A*)] and S — sin[so ln(A/A^)]. 
In the renormalization prescription for 173 in Eq. (j9bp . H 
must be a log-periodic function of A. The analytic ap- 
proximation for H derived in Ref. 16| is Eq. (|9c| with 
ho = 1. We find that to within our numerical accuracy 
of about 10^'', H is given by this analytic expression 
multiplied by the numerical constant ho = 0.879. The 
renormalization scale A* introduced by this prescription 
differs from by a multiplicative factor that is only 
known numerically: soln(A*/«:*) = 0.971 mod tt 11|. 
The function J in Eq. (|9b[) is not essential for renormal- 
ization, but it is needed to derive the Energy Relation in 
Eq. Q. 

Using the operational definition of C'2 and C3 in 
Eqs. H]) together with the Feynman-Hellman theorem, 
we can identify the 2-body and 3-body contact densities 
in the quantum field theory: 
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where H' and J' are the derivatives of H and J with 
respect to ln(A/A*). The contacts C2 and C3 are ob- 
tained by integrating these densities over all space. We 
have used the identity {ad/da)g2 — g'^f&'n'a as well as 
the expression for 53 in Eq. (I9b|) . The condition that 
C2 and C3 have finite limits as A 00 implies that the 
matrix elements (d^d) and (tH) scale as A^ and A^, re- 
spectively. Thus the (tH) term in Eq. pOap and the J' 
term in Eq. (jlOb[) can be omitted unless they are multi- 
plied by a factor of A. 

To derive the Energy Relation in Eq. (U), we express 
the interaction energy U — Jd'^R {Hint) in terms of the 
contacts defined by Eqs. ([TT 



(11) 
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The coefhcient of C2 scales as A as A — > cx), while the 
coefficient of C3 is a log-period function of A. The sub- 
traction terms proportional to C2 and C3 in the momen- 
tum integral on the right side of Eq. (jH) can be evaluated 



explicitly in terms of the ultraviolet cutoff A. After sub- 
tracting these terms from U , we find that the dependence 
on the ultraviolet cutoff A cancels. The subtracted ex- 
pression for U reduces to the last two terms proportional 
to C2 and C3 on the right side of Eq. This proves 
the Energy Relation and determines the constants A, ip, 
and fo in terms of the coefficients ho, jo, ji, and j2- 

The tail of the momentum distribution in Eq. ([2]) can 
also be derived by using the short-distance operator prod- 
uct expansion (OPE) |15| . The momentum distribution 
at wavevector k can be expressed in the quantum field 
theory as 



n{k)=Jd'Rjd' 



-ik-r 



{^\R~\r) V(H+ir)). (12) 



The coefficients in the OPE for ^li'^ and ■0 at equal times 
can be determined by matching matrix elements between 



asymptotic few-body states [15|. Alternatively they can 



be determined by matching Green functions in the few- 
body sector. The simplest choice of Green functions for 
the matching are those that are one-particle-irreducible 
(IPI) with respect to the atom field -0 and the diatom 
field d = The resulting OPE at large wavevector k 
can be expressed as 



/dV e-'''-'' V^(H-ir) ipiR+^r) 



fc4 4 



(13) 



where F{k) is the log-periodic function in Eq. Q and the 
additional terms are all suppressed by at least k^^. The 
Wilson coefficient of c?^c? was determined analytically by 
matching the diatom Green function. The Wilson coeffi- 
cient of tH was subsequently detemined by matching the 
atom-|-diatom Green function. This Green function can 
be expressed as a sum of loop diagrams involving the con- 
nected atom-|-diatom Green function, which can be calcu- 
lated by solving the Skorniakov-ter-Martyrosian integral 
equation numerically. Inserting the OPE in Eq. ()13p into 
the expression for the momentum distribution n{k) in 
Eq. ([T^ , we obtain the result for the tail of the momen- 
tum distribution in Eq. ([2]). Our direct calculation gives 
constants A and (j) that agree to within a few percent 
with the precise results given after Eq. ([3|). 

In a quantum field theory, rf transitions of an atom to 
a different spin state can be represented by an operator 
'i/'JV'i where i/^l creates an atom in the second spin state. 
The rf transition rate can be expressed as 

r{uj) = Imz/di e'("+")7d3i?/dV 

x(TV't02(i?+ir,O 4^{R-^r,0)). (14) 

We assume for simplicity that the scattering length for 
the second spin state and the pair scattering length for 
the first and second spin states are negligible compared 
to a. We can therefore take 1/^2 to be a noninteracting 
field. The rf transition rate at large w can be determined 
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by using the short-time operator product expansion for 
'0^''02 and -tplip Q. The leading terms in the OPE at large 
w can be expressed as 



+ 
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where -Frf('^) is a dimensionless function and the addi- 
tional terms are all suppressed by at least w^^/^. The 
Wilson coefficients for f/jti/; and d^d were determined an- 
alytically by matching the atom and diatom Green func- 
tions, respectively. The Wilson coefficient of tH was 
subsequently determined by matching the atom+diatom 
Green function. The function F^f{ijj) in Eq. (fT5|) has the 
form 

F,f{uj) ^ Do + Di\n{{-mojf'^/K) 

-fi:'2sin2[soln((-mcj)i/VK*) + 0rf], (16) 

where A is the ultraviolet momentum cutoff. Note that 
-Fi.f(a;) is logarithmically ultraviolet divergent. The nu- 
merical constants in Eq. ((T^ are Dq = 0.670, Di = 
-2.94, D2 = 1.16, and 0i.f = 1-33. To obtain the high- 
frequency tail in the rf transition rate, we insert the OPE 
in Eq. (fT5|) into the expression for r(a;) in Eq. (fT4|) and 
take the imaginary part. Terms that are analytic func- 
tions of u), such as 1/uj and 1/w^, do not contribute to 
the imaginary part. The leading terms at large uj are 
given by Eq. where G'rf(a;) — 2 ImFi.f(a; + ie). Using 
our result for F^f{u}) in Eq. p^ . we obtain the result for 
G,i{Lo) in Eq. 0. 

There are universal relations involving the 3-body con- 
tact for any system in which the Efimov effect arises in 
the 3-body problem. Another such system consists of 
fermions with three spin states, which we label 1, 2, and 
3. In the zero-range limit, the interactions are described 
by three large pair scattering lengths 012, ai3, and 023 
and an Efimov parameter . The discrete scaling factor 
has the same value e^^^° « 22.7 as for identical bosons. 
There are 2-body contacts C12, C13, and C23 associated 
with each pair of spin states and a 3-body contact C123. 
The tail of the number distribution 7ii(fc) for spin state 1 
is given by Eq. ^ with C2 and C3 replaced by C12 -I- C13 
and C123. An operational definition of C12 is given by 
Eq. (jlap with a and C2 replaced by ai2 and 2Ci2. An 
operational definition of C123 is given by Eq. (Ilbp with 
C3 replaced by C123. It is straightforward to derive the 
analogs of the Energy Relation in Eq. (|4|) and the uni- 
versal relation for the rf transition rate in Eq. (jS]). 

Many-body systems consisting of identical bosonic 
atoms or of fermonic atoms with three or more spin 
states are unstable due to recombination into deeply- 
bound dimers. The rates for these loss processes scale 



as for large a, which makes it difficult to test uni- 
versal relations for global equilibrium properties of the 
system, such as the Virial Theorem in Eq. ([5]). However, 
universal relations that govern the short-time behavior 
of the system, such as the tails of the momentum dis- 
tribution in Eq. ^ and the tail of the rf transition rate 
in Eq. ([6]), can be tested experimentally by using short- 
time probes of ultracold atoms, such as those that have 
already been applied to fermions with two spin states 
[lo| . These universal relations involve log-periodic func- 
tions, so they provide a new probe of Efimov physics in 
ultracold atoms. 

This research was supported in part by a joint grant 
from the Army Research Office and the Air Force Office of 
Scientific Research and by a grant from the Department 
of Energy. 



* Electronic address: 'braaten@mps.ohio-state. edu| 
^ Electronic address: kang@mps.ohio-state.edul 

* Electronic address: platter@chalmers.se 

[1] S. Giorgini, L. Pitaevskii, and S. Stringari, Rev. Mod. 
Phys. 80, 1215 (2008); W. Ketterle and M.W. Zwier- 
lein. Making, probing and understanding ultracold Fermi 
gases, in Ultracold Fermi Gases, ed. M. Inguscio, W. Ket- 
terle, and C. Salomon, (lOS Press, Amsterdam) 2008. 

[2] E. Braaten, arXiv: 1008.2922. 

[3] S. Tan, Annals of Physics 323, 2952 (2008); ibid. 323, 

2971 (2008); ibid. 323, 2987 (2008). 
[4] M. Punk and W. Zwerger, Phys. Rev. Lett. 99, 170404 

(2007) ; G. Baym, C.J. Pethick, Z. Yu, and M.W. Zwier- 
lein, Phys. Rev. Lett. 99, 190407 (2007). 

[5] W. Schneider and M. Randeria, Phys. Rev. A 81, 

021601(R) (2010). 
[6] E. Braaten, D. Kang and L. Platter, Phys. Rev. Lett. 

104, 223004 (2010). 
[7] F. Werner, L. TarrueU, and Y. Gastin, Eur. Phys. J. B 

68, 401 (2009); S. Zhang and A.J. Leggett, Phys. Rev. A 

79, 023601 (2009). 
[8] D.T. Son and E.G. Thompson, Phys. Rev. A 81, 063634 

(2010); H. Hu, X.-J. Liu, and P.D. Drummond, Eu- 

rophys. Lett. 91, 20005 (2010); W.D. Goldberger and 

LZ. Rothstein, arXiv: 1012.5975. 
[9] E. Taylor and M. Randeria, Phys. Rev. A 81, 053610 

(2010); T. Enss, R. Haussmann, and W. Zwerger, Ann. 

Phys. (N.Y.) 326, 770 (2011). 
[10] E.D. Kuhnle et al., Phys. Rev. Lett. 105, 070402 (2010); 

J.T. Stewart et al., Phys. Rev. Lett. 104, 235301 (2010). 
[11] E. Braaten and H. W. Hammer, Phys. Rept. 428, 259 

(2006). 

[12] F. Werner and Y. Gastin, arXiv: 1001 . 0774. 

[13] R. Gombescot, F. Alzetto, and X. Leyronas, Phys. Rev. 

A 79, 053640 (2009). 
[14] F. Werner, Phys. Rev. A 78, 025601 (2008). 
[15] E. Braaten and L. Platter, Phys. Rev. Lett. 100, 205301 

(2008) . 

[16] P. F. Bedaque, H. W. Hammer and U. van Kolck, Phys. 
Rev. Lett. 82, 463 (1999). 



